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We give analytical and numerical solutions for the zero modes of the Dirac operator in topological 
SU(2) gauge backgrounds at nonzero chemical potential. Continuation from imaginary to real 
chemical potential is used to systematically derive analytical zero modes for calorons at arbitrary 
holonomy and, in particular limits, for instantons and dyons (magnetic monopoles). For the latter 
a spherical ansatz is explored as well. All the zero mode exhibit stronger peaks at the core and 
negative regions in their densities. We discuss the structure of the corresponding overlap matrix 
elements. For discretized calorons on the lattice we consider the staggered operator and show that 
it does possess (quartet quasi-)zero modes, whose eigenmodes agree very well with the continuum 
profiles, also for nonzero chemical potential. 



I. INTRODUCTION 

The thermodynamics of Quantum Chromodynamics 
(QCD) has attracted a lot of attention in recent years, 
both experimentally and theoretically. While at high 
temperature the transition to a phase of deconfined glu- 
ons and quarks obeying chiral symmetry is rather well un- 
derstood, the exploration of a finite quark density is still 
immature. Theoretical progress in this part of the QCD 
phase diagram has mainly been hindered by the failure 
of the most powerful nonperturbative method from first 
principals, lattice gauge theory. At finite chemical po- 
tential the Dirac operator is not purely antihermitian 
and its complex determinant prevents an interpretation 
as a probability weight. This is the famous sign-problem. 
Ways to proceed include reweighting and expansion tech- 
niques on the lattice, complex Langevin dynamics, and 
studies of finite density in diagrammatic and related ap- 
proaches and approximations to QCD or similar theories 
such as random matrix theory, for reviews see e.g. [IIH]. 

Here we investigate the semiclassical approach to 
QCD, which is based on solutions of the Euclidean Yang- 
Mills equations. At zero temperature these are the cel- 
ebrated instantons [3^ , localized lumps of action density 
characterized by an integer topological charge. At fi- 
nite temperature, the corresponding solutions are called 
calorons [4]. Calorons have been found to be governed 
by an additional parameter, the holonomy, and consist of 
constituent dyons (magnetic monopoles, see below) [5l[6]. 
Individual dyons constitute static solutions at finite tem- 
perature. All these configurations are selfdual or anti- 
selfdual and thus minima of the Yang-Mills action in the 
given topological sector. 

The semiclassical theory of instantons and calorons 
with trivial holonomy has been studied at zero and finite 
temperature [7HTT] as well as nonzero chemical potential 
[HHH] (see [12 for a review) . The corresponding mecha- 
nism of chiral symmetry breaking is rather robust: exact 
fermionic zero modes become low-lying modes and thus 
generate a nonzero density of states at zero Dirac eigen- 
value, which via the Banks-Casher relation is linked to 
the chiral condensate. Confinement, on the other hand, 
could not be obtained from instantons nor calorons of 



trivial holonomy. This changes drastically for calorons 
at non-trivial holonomy, as has first been shown semi- 
analytically in [THj. The lesson is that the holonomy 
plays an important role with the confinement mecha- 
nism becoming close to the old Dual superconductor pic- 
ture. Promising results in this context are the possibility 
to show confinement analytically from the moduli space 
metric of dyons and even from a non-interacting gas 
of dyons [20^ . Although the interactions between calorons 
or dyons are not fully worked out, there are more hints on 
confinement, namely from the quantum weight I21|, the 
detection of calorons in lattice quantum configurations 
P^HMj and the vortex content of calorons [2S]. In the 
high temperature phase either dyons of particular kinds 
are suppressed 157] or dyons form molecules due to 
their fermionic interactions |28j . 

Instantons, calorons and dyons also share the support 
for fermionic zero modes, which are present due to the in- 
dex theorems . These modes localize to the back- 
ground as can be seen from the existing explicit solutions 
|33H35j . Generalizations of the index theorem to nonzero 
chemical potentials have been discussed in jSSl [37] . Our 
results confirm these as far as instantons/calorons of unit 
charge and dyons of fractional charge are concerned. For 
the semiclassical approach at finite density it is necessary 
to know the zero eigenmodes, both for the formation of 
the condensate and the fermionic interaction of the ob- 
jects. 

In this paper we present analytical and numerical re- 
sults for SU(2) zero modes at nonzero chemical potential. 
We show that the analytic continuation from imaginary 
chemical potential provides exactly the particular (bi-) 
orthonormalization necessary for the eigenmodes at real 
chemical potential. Together with the results from the 
ADHM-Nahm formalism |35| we systematically ob- 
tain the zero modes in caloron backgrounds at arbitrary 
holonomy and chemical potential. In particular limits we 
confirm the zero modes of trivial holonomy calorons and 
instantons known in the literature [^DHl^ as well as for 
dyons, in which case we show that the spherically sym- 
metric ansatz [28, 33J works. The typical features of these 
modes are stronger peaks at and oscillations away from 
the core, with the period of oscillation c>c cf. Fig. IT] 




FIG. 1. Zero mode in a background of a caloron with maximally nontrivial holonomy and /i — AT (colored, magnitude rescaled 
to its third root) together with the caloron action density revealing two constituent dyons (grey, cut along the caloron axis). 



This will have a wash-out effect in the fermionic interac- 
tions of these objects, which we briefly discuss. 

We also present numerical results for calorons dis- 
cretized on the lattice utilizing the staggered Dirac op- 
erator. Although the latter has four tastes and no exact 
zero modes of definite chirality, we identify modes smaller 
than the rest of the spectrum by several orders of magni- 
tude. These quasi-zero modes persist at nonzero chemi- 
cal potential, too. Taking the taste quartets together, the 
profiles agree very well with the profiles of the continuum 
zero modes. 

This paper is organized as follows: In the next sec- 
tion we discuss properties of the Dirac operator and its 
eigcnmodes at nonzero chemical potential. After that in 
Sec. |III|we perform the analytic continuation from imag- 



inary to real chemical potential. Sec. IV contains the 
corresponding analytic results for zero modes in dyon, 
caloron and instanton backgrounds. Sec.|V]is devoted to 
the lattice results for the staggered eigenmodes. Sec. |VI| 
gives first results for the overlap matrix elements of dyon 
zero modes. We summarize our findings in Sec. |VII[ 



II. 



DIRAC OPERATOR WITH CHEMICAL 
POTENTIAL 



In the continuum Dirac operator the chemical potential 
is added to the covariant derivative 



(2.1) 



(for our lattice Dirac operator see Sec. VB). We will al- 
ways focus on the massless case and analyse eigenmodes 



lp{^i)'ipriix; ij) = A„(/i)?/>„(a;; n) , 



(2.2) 



in particular zero modes, Ao = (arguments will be ne- 
glected whenever convenient). 

At finite temperature T = 1//3, physical modes are 
antiperiodic in the temporal direction 

^(t-f/3,f) = -^(t,f), (2.3) 

but we will also consider modes periodic up to a phase. 



A. Symmetries 

Since the Dirac operator is chiral at arbitrary chemical 
potential, {-0(/.t),75} — 0, its eigenvalues always come in 
pairs ±A. 

While the Dirac operator at vanishing /i is antihermi- 
tian with its eigenvalues purely imaginary, this important 
property is lost for real However, for the hermitian 
conjugate operator the following relation holds, 

0^(m) = -0(-a^), A*(m) = -A(-a.), (Ai€R)(2.4) 

(for complex /i's replace — /i by —fi*). 

For gauge group SU(2), due to the pseudoreal nature 
of this group, the Dirac operator with (real) chemical 
potential obeys an anti-unitary symmetry 

lp*{^i)^-W^lp{fi)W {SU{2),fieR) (2.5) 

where W is unitary and /i-independenl[^ with W*W — 1. 
In Random matrix theory this symmetry qualifies iIp{Q) 
as belonging to the universality class of the Gaussian or- 
thogonal ensemble. Moreover, in a particular basis ilp{0) 
is real symmetric with real eigenmodes (see below; for the 
different anti-unitary symmetry of our lattice Dirac op- 
erator see Sec. VB|. In the next subsection, we will use 



this symmetry to show that the SU(2) zero modes have 
real densities. 

For operators without definite hermiticity one has to 
distinguish eigenvalues with corresponding left and right 
eigenmodes from singular values as will be discussed in 
the next subsections. 



^ In the Weyl representation one has W = 7270 ''"2, with the 
charge conjugation 7270 relating gamma matrices to their com- 
plex conjugates and the second Pauli matrix T2 relating the Pauli 
matrices as generators in color space to their complex conjugates. 
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B. Left and right eigenmodes 

The right eigenmodes of the Dirac operator are given 
through Eq. (2.2 1, whereas the left eigenmodes are de- 



fined analogously 

xL(Ai)0(M) = xL(M)Am(Ai) ■ 



(2.6) 



For our specific case, the left eigenmodes are related to 
the right eigenmodes at opposite fi 



(2.7) 



(which follows from the hermitian conjugate of Eq. (2.2) 



and the property (2.4)). Thus, one may replace one by 
the other. 

From the difference of the scalars product of xLIm) 
with (2.2) and of (2.6) with ipnifJ-) one immediately con- 
cludes that right and left eigenvalues agree and that 
left and right eigenmodes are orthogonal. Defining the 
(pseudo- ) density 

Pmn{x; m) = i^lnix; -/l)V'n(x; /i) 

= ^V™(a^;-/i)*C(a^;M) (2.8) 



where the index a stands for both color and spin, the 
x-integral over p„i„ is 



(2.9) 



Such a basis is called bi-orthonormal. At vanishing chem- 
ical potential this of course reduces to the orthonormal 
basis of eigenmodes of the purely antihermitian Dirac op- 
erator. 

Note that for /i 7^ the densities Pmm{x) of individual 
eigenvalues need not be real everywhere in space-time, 
although their integral is positive. 

Densities of SU( 2) z ero modes can be made real thanks 
to the symmetry (2.5). For this gauge group «0(/x) can 
be transformed to 5* + ^.A with real symmetric S and 
real antisymmetric A |43j . For a real matrix the zero 
modes can be chosen real (as the complex conjugate of 
a zero mode is a zero mode again). Hence the density 
Pm{x) is real in this and every other basis (as this fermion 
bilinear is an invariant). How the densities come out real 
in an arbitrary basis is shown in App. [Aj Still, negative 
densities can, and will, occur. 



C. Singular value decomposition 

Arbitrary matrices can be written in a singular value 
decomposition, for the Dirac operator 



[/t0F = diag(ei,6,---), en>0 



(2.10) 



where U and V are unitary transformations and the en- 
tries ^„ - which are non-negative - are called singular 



values. The columns u„ of U and u„ of V are called left 
and right singular modes (vectors), respectively. They 
are orthonormal among themselves, but not necessarily 
to each other. 

For generic configurations at nonzero /x, these singular 
values and modes are not related to the eigenvalues and 
eigenmodes of the operator Ip itself, rather to those of 
its product with 0^: the left and right singular modes 
are the eigenmodes of the hermitian operators Iplp^ and 
0^0, respectively, and the nonzero singular values ^„ 
are the square roots of the eigenvalues of these operators 
(which coincide between Iplp^ and Ip^ Ip). 

Zero singular and eigenvalues, however, constitute an 
exception as the zero modes are identical, vq = "00 £^nd 
""0 = Xo Thus, in order to determine zero modes of 
the Dirac operator, both the eigenvalue and the singular 
value method can be used. 

At the fundamental level, the singular value decompo- 
sition is slightly more powerful, as it always exists and its 
complete and orthonormal basis can be used to decom- 
pose the fermion fields in the path integral. Eigenmodes, 
on the other hand, do not always provide a complete ba- 
sis. The simplest example for a matrix of this kind is 

, a 4 X 4 example is discussed in [37'. Finite /i 

Dirac operators with this property occur in principle, but 
the corresponding configurations are non-generic. 



D. Index theorem 

We start by recalling the case of vanishing In a 
chiral basis the Dirac operator can be written as 




(2.11) 



Then 75 — diag(l2,— 12) and left and right handed zero 
modes consist only of lower and upper components, re- 
spectively. Let us denote their numbers by tll^r — 
dim ker Dj^ and specialize to positive topological charge 
Q for simplicity. 

The index theorem reads ri/j ^ — Q. It is ful- 
filled by generic configurations through possessing Q 
right handed, but no left handed zero modes. 

With nonzero chemical potential the index theorem has 
been discussed in Refs. j36j|37j. It now also incorporates 



^ That a right zero cigenmodc of is also a right zero eigenmode 
of and thus a zero singular mode follows immediately; the 
opposite direction of the proof uses that 0^0'(/> = induces 
l^i/ip = 0, which means that Ipip must vanish. The same argu- 
ment applies to left zero eigenmodes. 
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the hermitean conjugates of D^ ji, for which the relatfon 
(2.4) impHes 



Dl{^l)^-DRi-^l), = -DiHi) (2.12) 

A straightforward generaUzation of the index theorem to 
real chemical potential has been obtained in |36) 

nB.i^J■) - nL{^J■) = Q , nii(-n) - nL{-^l) = Q (2.13) 

where we have trivially added the index theorem at op- 
posite ^. The more appropriate index theorem claimed 
in |37j (Eq. (D.3)) - derived via the singular value de- 
composition - translates via (2.12) into 



nR{-y)~nL{ii)^Q (2.14) 



The two index theorems agree for almost all configura- 
tions [33. 

Our background configurations obey both index theo- 
rems in the generic way by possessing one right handed, 
but no left handed zero mode. Due to the latter, they are 
not sensitive to the difference of the two index theorems. 



III. ANALYTIC CONTINUATION IN THE 
CHEMICAL POTENTIAL 

In this section we will demonstrate that eigenmodes at 
imaginary chemical potential can be analytically contin- 
ued to yield eigenmodes at real chemical potential obey- 
ing the correct bi-orthonormalization. 

Let us consider the operator 

Ip'{^) = -1,D, + i-io^T = Ip{-i^T) (3.1) 

which is nothing but the Dirac operator with an imag- 
inary chemical potential —i(pT. The minus sign in the 
argument on the right is for later convenience. At real 
(f the operator Ip is purely anti-hermitean, thus in 
its eigenequation the eigenvalues are purely imaginary 
and the eigenmodes are orthonormal in the conventional 
sense, i.e. the densities 



cf. Eq. ( |3.1[ ). Then, analytically continuing the 
eigenequation, the ip'„{x;i^/T) obey the eigenequation 
of 0(^) and are thus candidates for eigenmodes. 



(3.5) 



That these modes are indeed eigenmodes, all that needs 
to be checked is their bi-orthonormalization. For that we 
carefully perform the analytic continuation in the densi- 
ties 

a 

= EC(-Wr)*C(*Ai/T) 

a 

= J2i^'^(.-^^rrni^^)^Pmni^i) (3.6) 

a 

valid for all x (suppressed). The opposite /i-arguments 
of the two factor^ are exactly what is needed to match 
the pseudo-density pmn (m) at real chemical potential, cf. 

Then the desired bi-orthonormalization of it, 
is inherited from the analytic continuation of 



Eq, 



Eq. (2.9 



2.i 



the conventional orthonormalization (3.3 1 of eigenmodes 
of the purely antihermitian operator . Moreover, we 
have started with eigenmodes ip' being antiperiodic in xq 
and so are their analytic continuations; actually all kinds 
of boundary conditons are continued to the eigenmodes 
at real chemical potential. Hence, the modes ipni^) ana- 
lytically continued to (p = i/i/T are indeed the eigenmodes 
at real chemical potential, i.e. Eq. (3.5) does hold. This 



continuation applies to all spin and color components of 
the eigenmodes and, of course, also to the eigenvalues, 
i.e. A'„(z/i/r) = A„(/i). 



The dimensionless parameter Lp in Eq. (3.1 ) can be in 



terpreted as a constant trace part of the gauge field and 
can also be transfered to the temporal boundary condi- 
tion defining 



M^-^)=e^^'/^^'^{x-^) 



(3.7) 



which are eigenmodes of the Dirac operator at vanishing 
chemical potential, but periodic up to a phase (p (on top 
of the antiperiodicity) 



are non-negative at every space-time point x and inte- 
grate to 



d'^x p'^„{x;ip) = S„ 



(3.3) 



Now we analytically continue (p to imaginary values 
ip/T assuming this is possible continuously. The Dirac 
operator becomes that at real chemical potential. 



y{tf,/T) = 0(/i) 



(3.4) 



e "^■0„(i,x;.^) , ip = Tr-ip (3.8) 

where we have also defined the total phase ip for later 
use. The pseudo-densities of these mode still agree with 
those of the original eigenmodes, p = p'. 



^ In going from the first to the second fine in \3.(>\ we have used 
that tp' is a smooth function in real ip with Taylor expansion 
^'{(p) = Yln'^^'f'" ^'^"^ complex conjugate reads ^'*{<fi) = 
En '^n'?"- analytic continuation of the latter we plug in 

an imaginary if, ^" {v)\f,=i^/T = SZ„ < (^m/^")", which agrees 
with the complex conjugate of ip'{—ifi/T). 
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In particular, knowing the zero modes in an analytic 
fashion for (small) imaginary chemical potentials, one im- 
mediately gets the zero modes of the desired Dirac oper- 
ator at real chemical potential. In Sec. |IVC| we will make 
use of this analytic continuation, since in caloron back- 
grounds (and limits thereof) the zero modes are known 
analytically for whole intervals of imaginary chemical po- 
tentials/boundary conditions. 



IV. ANALYSIS OF ZERO MODES IN SPECIFIC 
TOPOLOGICAL BACKGROUNDS 

A. Single dyon zero mode 

Here we discuss the solution of the zero mode on a 
single dyon for arbitrary boundary condition and with 
arbitrary chemical potential /i. By dyon we mean a mag- 
netic monopole obeying the BPS condition, which in a 
3-|-l-dimensional language amounts to possessing an (Eu- 
clidean) electric charge of same magnitude]^ The deriva- 
tion of the zero mode makes use of the spherical symme- 
try and closely follows [IHl 123] ■ 

We work in the BPS limit and in a hedgehog gauge. 
The (anti)self-duality equations for the gauge fields are 
solved for by (see e.g. [IS]) 



4" = Aea^kr'' , A 



^0 



U 



1 - 



VT 



sinh(wr) 
vr coth(wr) 



(4.1a) 



±i;, (4.1b) 



with f = x/r being the unit vector in the radial direction 
and = ^"iV "^ith t° the Pauli matrices. The up- 
per sign refers to the sclfdual and the lower to the anti- 
selfdual field equationtj^ and v is the asymptotic value of 

^0, 



- lim J{A-) 



(4.2) 



playing the role of a "Higgs" vacuum expectation value. 
For simplicity we will take u > 0. The vev also enters 
the action and topological charge of the dyon, which we 
define in the four-dimensional language: 



S 



1 



/3 



Att/Sv , Q 



S 

8^ 



2ttT 



(4.3) 



It was argued in |44j that the term "dyon" is not entirely appro- 
priate for these objects. They suggest that the term "monopole- 
instanton" is more accurate. Nevertheless we will use the term 
dyon in this work to mean the constituent of a finite temperature 
instanton, see below. 

We define selfduality as _Fpj/ = f-fivprrFva l"^, where iJ,,v,p,a = 
0, ... ,3. From eoi23 = 1 this yields, e.g. Fqi = -F23. Note that 
self-duality equations with fi,u, p,cr = 1, . . . , 4 and £1234 = 1 
yield F41 = -F23. 



We restrict ourselves to u < 27rT, such that Q G [0, 1]. 

For the zero modes we make the usual radial ansatz 
for Dirac spinors (see ^51 1551 HS] ) 



(4.4) 



where A is the color index and a the spin index, and 
e = ?T2 is a completely antisymmetric rank-2 tensor. Ob- 
viously these functions are periodic up to a phase e"*"^, 
cf. (3.8). We stress that the gauge background is a - 
so far - i-independent embedding of a three-dimensional 
dyon into x R-^. 



The Dirac equation (2.11 for left and right-handed 
spinors is 

iicTnMD^)AB - ^^SABSc.p){i^R)p = , (4.5a) 
{{anMD^)AB - MAB<5o0)(V'L)f = . (4.5b) 

Here cr^ = (1, — «t*), ct'^ = (1, ir*) and we get the follow- 
ing equations for the functions ai^2{i~) 



— H 

dr 2 

da2ir) 



ai — + iipT)a2 = , (4.6a) 



f±H-2A 2\ . rr., 

, ^ + -ja2-i{fi + i^T)ai = . 

(4.6b) 

where the upper sign is to be taken for the left, and the 
lower for the right handed spinor component. One can 
see once more that the boundary phase is equivalent 
to an imaginary chemical potential iipT. In the following 
we will therefore use 



= ^ + iipT 



(4.7) 



To discuss the normalizability of solutions, let us spe- 
cialize to = for the moment. While A decays for 
large distances, % approaches v for self-dual and —v 
for anti-self-dual dyons, respectively, see Eq. (|4.1b|). In 



the asymptotic linear system of differential equations 
(4.6) the chemical potential gives rise to trigonomet- 
ric functions. The "H-term leads to exponential func- 
tions and thus determines, whether the solutions are 
normalizable. For the latter the equations must read 
da,{r)/dr + vaj2 + . . . = 0. Hence a normalizable solu- 
tion for the left and right spinor exists if the background 
field is self-dual or anti-self-dual, respectively, and the 
solutions behave like e^™'!'^ asymptotically. 
The full equations are solved by [28] 



ai,2(?';M) = c 



^ vr sinh(ur) 



Xi[fv) — 2 — sin(r/i^) + tanh(wr/2) cos(r/i^) 



(4.8a) 



(4.8b) 



X2{rv) = =F« 2 — cos{rfi^) - coth(i;r/2) sin(r/i^) . 

(4.8c) 
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FIG. 2. Densities poo{r', f^) of periodic dyon zero modes, 



Eq. (4.81 with (y3 = 0, as a function of rv for /i = 0.5v (solid 
blue), jj, — V (dashed red) and /i — 1.2v (dot-dashed purple). 



FIG. 3. Logarithmic plot of the absolute value of the density 
Poo(r; jj.) for fj, = 3v (solid) and that of the envelope Eq. (4.13 1 
(dashed). 



In these formulae, the boundary condition if has been 
reintroduced, as the imaginary part of ^^p, cf. (4.7). It 
leads to exponential parts e''''^'', which compete with the 
exponential decay e""''/^. Consequently, the solution is 
normalizable only if ip £ {—v/2T,v/2T), i.e. around pe- 
riodic boundary conditions. Later we will consider mod- 
ifications leading to antiperiodic zero modes. 

Since we are interested in the dependence on ^ we will 
set again ip — in what follows, i.e. analyze periodic zero 
modes. 

As we discussed in Sec. |II B[ the relevant density of 
these modes is given by = poo(p), 

Pooir^n) ^ al(r; -fi)ai{r; fi) + a2{r; -fj,)a2{r; IJ.) (4.9) 



The normalization in the sense of J d'^x poo = 1, 



cf. 



Eq. (2.9), fixes the constant 



(4.10) 



From Eq. ( |4.9[ ) one can already see the possibility for 
this density to become negative, through a2- Consis- 
tently, at vanishing fj, the (/i-odd) function a2 vanishes, 
and the density is positive definite. Profiles of periodic 
zero modes for several values of the chemical potential 
are plotted in Fig. [2j 

The solutions clearly have an exponential envelope, of 
the asymptotic shape e~^^^'^ /^/r, see Fig. ^ For the 



prefactor in the envelope one has to consider |Xi T ~ 1X2 1 
For large vr one can set tanh(vr/2) = 1 = coth(wr/2) up 
to exponentially small corrections and get 



\Xi{rv)\'' ~ IX2MP = ( 1 + 4 ^ ) sin{^{^^/v) + 2H 



S{ii/v) = arctan 



1 - Afi^v^ 
Aijl/v 



(4.11) 
(4.12) 



The factor in the parenthesis bounds this term and from 
Eqs. (IXsl |4.10|) we get 



(4.13) 



again up to exponential corrections in vr. This envelope 
has been included in Fig. [3| 

Another property of the zero mode is that its value at 
the dyon core grows with the chemical potential, which 
can also be seen in the examples of Fig. [2] Taking the 
limit r in the solution Eq. (4.8) it is easy to see that 



Poo(0;m) 



T 
IGttw 



(4.14) 



Probably the most characteristic feature of the zero 
mode density are its negative regions. For any nonzero 
chemical potential there are infinitely many negative re- 
gions (recall that the far away regions have exponentially 
small densities though) , separated by infinitely many ze- 
roes. This is very obvious in a logarithmic plot like Fig. |3] 

For /X <C w we can find these zeros explicitly by de- 



manding that Eq. (4.11) vanishes 



1 

'2p 



5{plv) 



Tr{n + 1) 
24J. 



TT 



2fi ' 



(4.15) 



where we expanded S{ii/v) for small fi. For large fi we 
expect the first zero to be small, i.e. ru <C 1, and get 
from Eq. (IXsl) that 



IXil 



IX2I 



-2 + (2 - AfPr'^) cos(2^r) + 4^r sin(2/ir) 

(rv)^ 

(4.16) 

which has a zero for ro/x = 2.04279. Fig. |4] shows the 
numerical solution of the first zero, as well as the limits 
of p/v 1 and fi/v <C 1. 

Finally, we discuss the possibility of zero modes with 
physical, i.e. antiperiodic, boundary conditions. The 
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FIG. 4. The first zero ro of the density poo(''; /i), as a function 
of (in units of v) ~ black line. The dashed line and the dot- 
dashed line show the first zero for small fi given by Eq. ( |4.15[ ) 
with n = and for large fi given by ro — 2.04279/^, from 



Eq. (4.161, respectively. 



static dyon discussed so far has normalizable zero modes 
for near periodic boundary phases, namely 



V V \ 

e ( - — , —j 



(4.17) 



There is another embedding of the dyon into x R^, 
that will play a role in the caloron solutions as well (as the 
complementary constituent dyon). To see how it arises 
note that the solution Eq. ( 4.8 1 has the topological charge 



Q = v/{2TrT), which is fractional for v/T G [0,27r]. In 
order to get a solution which complements the topological 
charge, we should replace the parameter v v — 2TrT—v 
and obtain a topological charge Q = 1 — Q. However, this 
change of parameters alone is not sufficient to match the 
two configurations. For 3 + 1-dimensional configurations 
it is useful to define the asymptotic Polyakov loop, also 
named holonomy. 



( ^\ 

Pqo = lim "Pexp i / dtAQ{t,x) 

\ Jo J 



(4.18) 



with V denoting path ordering along xq. The starting 
configuration with gauge fields given by Eq. ( |4.1[ ) has 

itrPoo =cos(v/2T) (4.19) 



Poo = exp (i— f • r 



whereas the Polyakov loop of the new solution would be 
P„o = exp {^^r --r) = - exp (^i^ 



exp^^— r-f) = -Poo 



and thus different. 

Therefore, consider the gauge transformatiorj^ 

t' 



ri(t, x) — exp jTrf • f 



(4.20) 



(4.21) 



It is antiperiodic, f2(/3,x) = —1 = — f2(0,x). As the 
gauge fields transform in the adjoint representation, they 
remain periodic. Hence, the dyon gauge field of ( |4.1[ ) with 
the replacement v ^ v transformed with is a valid con- 
figuration on S^xR^. The gauge transformation does not 
alter its action and topological charge, so they are still 
proportional to v (cf. (4.3), both are gauge invariant), 
and thus also inherits (anti)selfduality; its gauge fields 
are just not static anymore. However, n does affect the 
Polyakov loop as Poo ~Poo = Poo- This means that 



the static dyon solution Eq. (4.11 and the dyon solution 
with V ^ V — 2-kT — V and the gauge twist of Eq. (4.21 ) 
both have the same asymptotic Polyakov loop Poo and 
complementary topological charges Q,Q = 1~ Q respec- 
tively. 

In calorons of integer topological charges, these dyon 
solutions occur together (see below) and are sometimes 
refered to as "light" and "heavy" dyon. The reason for 
this nomenclature is because the Polyakov loop branch in 
full QCD, i.e. when fermions are dynamical, interpolates 
between unity and (near) zero from high to low temper- 
atures. This would correspond to the parameter v in the 
range v £ [0,7rr]. In this range the untwisted dyon (i.e. 
the one without the gauge transformation Eq. (4.21 )) has 



the action and the topological charge proportional to v 
and the one with the gauge twist has the same propor- 
tional to V. Since in this range v < v the dyons untwisted 
and twisted dyons are referred to as "light" and "heavy" . 
Finally, we note that the fundamental fermions feel the 



gauge transformation of Eq. (4.21), they obtain an addi- 
tional factor of —1 in their periodicity property. Conse- 



quently, the parameter ip used in the ansatz Eq. (4.4 1 is 
not the boundary phase as defined in Eq. (4.21 ). Rather, 
we should replace ip ^ if — n — —Cp in this ansatz. With 
the additional replacement v ^ v \i follows that these 
zero modes on top of the twisted dyons are normalizable 
for 



pe 



2T' 2T/ ^ ^ V2T' '^^ 2t) ^ ' ' 



which is the complementary range of the static dyons, 
Eq. (4.171. For the chemical potential the replacement 
reads /z^ — > fi-tp — fi ~ iipT. Antiperiodic fermions are 
included as p — with /i remaining the same. 



B. Calorons and instantons 

In this section we will demonstrate that zero modes 
at imaginary chemical potential are included naturally 
in the ADHM-Nahm formalism [351 131113 , so that they 



In the gauge of the caloron, which combs the field Aq to the 



third color direction, the associated gauge transform is Q = 
exp{iiTT-j t/ P). For higher gauge group such gauge transforma- 
tions need to obey f2(i = /3) = zQ{t = 0), where z is an element 
of the center of the gauge group. 
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are available - at least in principle - for all selfdual so- 
lutions of four-dimensional Yang-Mills theory. After an- 
alytic continuation one thus obtains zero modes at real 
chemical potential. We will demonstrate this explicitly 
for SU(2) calorons and then in particular limits recover 
the zero modes of dyons from the previous section and of 
calorons and instantons from the literature. 



1. Calorons, especially in SU(2) 

Calorons are (anti-)selfdual Yang-Mills fields ('instan- 
tons') at finite temperature T, i.e. over the space-time 
manifold x B? where the Euclidean time is compact- 
ified on a circle with circumference 1/T. 

The asymptotic Polyakov loop plays an important role 
as a parameter of caloron solutions. For the gauge group 
SU{2) we parametrize it with the holonomy parameter 
we [0,1/2]: 



Poo = exp(27riwT3) , ^ tr -Pa 



cos(27rw) . (4.23) 



The two extreme cases uj = 0, uj — 1/2 amount to 'triv- 
ial holonomy' Poo = ±l2, while the middle w = 1/4 
with tr Poo = is refered to as 'maximally nontrivial 
holonomy'. From the corresponding asymptotics of the 
Polyakov loop, Eq. (4.19 ), one can read off the connection 



to the (dimensionful) "Higgs vev" of Eq. (4.2) 



V = AttojT, 



4ttT 



(4.24) 



As a new aspect in comparison to zero temperature 
instantons (which can be thought of as coming from di- 
mensional reduction and symmetry breaking by Poo), 
calorons in gauge group SU {Nc) with topological charge 
Q consist of A^dQI constituent dyons. These appear as 
lumps of topological density with fractional topological 
charges given by the holonomy. In SU (2) the topological 
charges are 2a; (indeed, the identification above yields the 
monopole topological charge of Eq. (4.3)) and the com- 
plementary 1 — 2aj = 2d). We will refer to the dyon with 
the higher topological charge as the "heavy dyon" and to 
the other as the "light dyon" . The positions of them are 
moduli of the solution. 

The explicit gauge field of SU{2) calorons have first 
been given in [51 [B] (for multi-calorons see |1S]). We 
follow the ADHM inspired formalism of [5], reinstating 
the temperature T. The distance to the dyon of masses 
uj and Lo will be denoted by r and s, respectively, and the 
distance between the dyons hy d = np^ /(3. In the limit 
of p — > CX3 a single dyon will be recovered. 

Note also, that the calorons in these references come 
in a particular gauge with two aspects. The first one, 
refered to as algebraic gauge, means that gauge fields are 
periodic up to the holonomy and — >■ far away from 
all the dyons. In other words, the holonomy Poo has 
been put into the twisted periodicity conditions. This is 
useful since such asymptotics simplify the superposition 



of dyons within the caloron and also the approximate 
superposition of dyons and antidyons in ensembles. The 
gauge transformation back to periodic gauge is given by 
(Poo) which is non-periodic. 

Secondly, the caloron is in a quasi-stringy gauge. For 
magnetic monopoles/dyons it is well-known that when 
one tries to "comb" the Higgs field Aq to asymptotically 
approach a constant color direction, a singularity along 
a Dirac string appears. This is so because the winding of 
Aq, which equals the magnetic charge and is apparent in 
the hedgehog gauge, constitutes a topological obstruction 
against such a stringy gauge. In the caloron solutions, Aq 
is proportional to T3 almost everywhere; the correspond- 
ing winding caused by each dyon is performed in an ex- 
ponentially small solid angle around the line connecting 
the dyons |25j, cf. Eq. ( |4.45 1 below. This quasi-stringy 
gauge leads to large gradients along this line, but no sin- 
gula rities (a part from the caloron's center of mass, see 
SeclIVDSIand SecM. 



C. Profiles of caloron zero modes 

The fermionic zero mode of the caloron (at = 0) 
hops between the constituent dyons as a function of the 
periodicity angle. This has first been demonstrated in 
j34j and is in accordance with the discussion about dyon 
zero modes in Sec. |IV A[ for near periodic zero modes 
the static dyon supports the zero mode, whereas for the 
complementary range of boundary conditions around an- 
tiperiodic the complementary time-dependent dyon sup- 
ports the zero mode. Index theorems on x R'^ can 
be found in |30H32| . This scenario will remain with real 
chemical potentials. The explicit formulas for the zero 
modes are included in the ADHM formalism and are 
given by the ADHM Greens function / (see below). We 
will analytically continue them to obtain zero modes at 
real chemical potential. 

In [34l [35] explicit expressions for the zero mode with 
periodicity up to a phase and the holonomy (stemming 
from the algebraic gauge) 



+ P,X)= e-2"-/^Poo*,/T(t, X) 



(4.25) 



were given (the dual variable to z comes with l/T since 
it has inverse length dimension). Gauging away Poo and 
comparing to (3.8) we identify 2ttz/T — tp, thus the zero 
mode with phase boundary conditions reads 



= (Poo)"*/^*l/2-W2.(2:) 



(4.26) 



The zero mode at imaginary chemical potential follows 
by virtue of (3.7) 



^P'oix;^) = (e^^Poo)-*/^^i/2-^/2.(a^) ■ 



(4.27) 



Finally, the zero mode at real chemical potential follows 
through analytic continuation, i.e. the replacement (f = 
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i^/T, which is the antiperiodic mode. Altogether, the results 

from [34l [35] have to be taken at a complex z = T/2 ~ 
i/j,/27r and multiplied by a real factor and color phases 
V'o(a;;/i) = e''*(Poo)^*/'^*i/2-W2^T(2:) (4.28) from P^o = exp(27riwT3). We obtain 



-2-Kiu) tT 



T 

i-^T, — — i — 
' 2 2n 



+ a./((i-.)r,|-.A 



1 0-1. 



with the following auxiliary functions 



1p = -Ct + CrCs + 
7-2 + + (p 



+ .S^ - 

2rs 



C* + CrCs + 



2rs 



(4.29) 

(4.30a) 
(4.30b) 



Ct =:cos(27rtT) , Cr = cosh(47rrwT) , Cs = cosh(47rsi:jr) , Sr = sinh(47rrwT) , Sg = sinh(47rsi:jr) (4.30c) 
/(z, z') = e^'^**^^^^' Vr(rs^^)-i{e-2''**'^ sign(^-^') ginh(27rr|z - z'|)s 

- cosh(27rr(z + z' - T)) dsCs + ^(s^ - + d2)s. 



+ cosh(27rr(2(I)r ~ \z ~ z'\)) dsc^ + ^{r^ + s'^ + 



+ sinh(27rr(2wr - \z - z'\)) 



(4.30d) 



We repeat that r and s are the distances of x to the dyon centers and d is the distance between the dyons. The caloron 
gauge field is given in terms of these functions, too The Greens function at the arguments needed simplifies to 



■'^ '2 27r' 



--IT itT iu>tT ^-Ht^rp 



2-iTitT 



sinh(7rrT — 2'KrLjT — ifir)s 



cosh(7rrT — 2TTrujT + ifir)rSs + sinh(7rrT — 2iTrujT + i/ir) 



ds. 



/((I - a.)T, I - *f ) = I - 

Z ZTT Z ZTT 



The last relation reflects the anti-unitary symmetry (2.5) of SU(2) and has also been used (at = 0) in 
the components of the zero mode. 

For the density a simpler formula applies 



1 



(2vr) 



2 1^ 



(4.31) 
(4.32) 
to relate 

(4.33) 



with the Greens function 



T . fi T . /i 

2 ^*2^' 2" ^*2^ 



^ = 7rT(rsV')"^| - r"^cos(2/^r) 



dsCs + -{s'^ - + d^)s. 



+ r 'Cr 



dsCs + -(r^ + + d^)s. 



sc., + dsc 



(4.34) 



The reader will note that the function /(z, z') given above and van Baal f5], we have given only the relevant part for 
is not the full function for arbitrary arguments z, z' . The our discussion valid for z/T, z' jT S [w, 1 — w]. 
full function can be found in the original papers by Kraan 

The resulting zero mode densities are shown in Fig. [5] 
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-2-10 1 2-2-10 1 2-2-10 1 2 



X//J X//S X//S 

FIG. 5. Densities of zero modes of a caloron (solid black line), single dyon (dashed blue), and single dyon with correction to 
the holonomy v ^ v + 1/d, where d is the distance to the other dyon (dot-dashed red). The plots are given for 3 values of 
/x/T = 0, 1, 5 from left to right. The densities are along the axis transverse to the position of the second dyon in the caloron. 
All plots are given for p = 0.6/3 and for maximally nontrivial holonomy lj — iL — 1/4 — v/{4nT) = v/{4:TtT). 



where we also compare with the single dyon zero mode 
discussed previously. Again these densities have nega- 
tive regions. Note an interesting effect that even for fairly 
large values of the size parameter p, i.e. for well-separated 
dyons, the two densities do not match very well. This 
mismatch can be compensated for by a trivial adjust- 
ment to the effective holonomy of a single dyon. Namely 
the solution for a single dyon with holonomy v should 
be changed to that with holonomy v + 1/d in order to 
match the solution in the caloron. In other words the 
holonomy at the dyon centre is not v, but is corrected by 
the asymptotic value of the other constituent of caloron 
and is approximately 11+1/ |d+r| fnv+l/d. This has been 
noticed in the full solution [3S] where it was shown that 
single dyon "mass" is renormalized by the presence of the 
other dyon. This renormalization of mass is not new and 
has been know for a long time for Yang-Mills+adjoint 
Higgs system and that it effects the classical interactions 
of monopoles [49H52] . Fig. |5] also shows the solutions 
with this adjustment. Considering the crudeness of this 
approximation, the agreement is remarkable up to fairly 
small instanton size. 



D. Various limits 

In this section we will discuss various limits of the gen- 
eral caloron solution of Eq. ( |4.29[ ) and ( |4.33|4.34| . In 
particular, we take the trivial holonomy limit and com- 
pare it to previous work, as well as the zero temperature 
limit and the large separation limit, which leads to single 
dyons. These limits are concisely summarized in Table |l] 



immediately leading to (writing out d — np^T) 



= 1 



-Ct + Cr , 



-Ct 



= n 



Tip'T 



Ct 



(4.36) 



where we have identified 4> with the function 11 used in 
the literature [13112 [S3]. 



Concerning the Greens functions Eq. (4.30d), all de- 



pendence on s drops as it should to restore spherical 
symmetry. From Eq. (4.34) / we obtain 



~ ttT ( d , /„ N 

J = — r 1 ( ^ cos(2^r) + Cr) + s 



)+Sr} 



, , , , , (4.37) 

rip L r J 

^2^2y2 _ cos(2^r) + cosh(27rrT) + ^ sinh(27rrT) 
~ r2 - cos(27rtT) + cosh(27rrr) -I- ^ sinh(27rrr) 

We recall that the zero mode density follows from act- 
ing with the four-dimensional Laplacian on this function, 
Eq. ([433|. 



For the zero mode components, the corresponding 
Greens function simplifies to 



T 



— f 

rip 

+ e 
2ttT 
rip 



-iTTtT 



sm 



h(7rrT + i/ir)} 



e '^'{cos(^r) cos(7rtT) sinh(7rrr) 



+ sin(^r ) sin(7rtr) cosh(7rrr) } (4.38) 



1. Trivial holonomy limit lu 

In the trivial holonomy limit uj — > 0, w — > 1/2 the 
auxiliary functions simplify in the following way 



Cr = cosh(27rrr) , Cg = 1 , 
Sr = sinh(27rrT) , = , 



(4.35) 



This agrees with /(T, T/2 — ip,/2Tr) because the Greens 
function f{z, z') is periodic in both arguments bycon- 
struction [S], and also because both a re real using (4.32 1. 
The projection matrices in Eq. (4.29) then add up to the 
unity matrix. Put together, the zero mode reads 

Mf^)^ = £e^'VlidJ (^0, 1 - {a,e)Ac. (4.39) 
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To compare to previous work we write 
2 Zir 11 



(4.40) 



$ = (n 1) ( cos(7rtr) sin(/ir)siii(7rtr) I 



sinh(7rrT) 



which turns the zero mode into the form given in 
for this case of trivial hofonomy. 



2. Instanton limit T — > 

In the zero temperature hmit, caforons become instan- 
tons irrespective of the original holonomy. Indeed, all lj 
and w enter with T such that they disappear in this limit. 
From the previous trivial holonomy case we get 



/ 



1 



r2 + <2 + p2 



1 



and 



n = 1 



$ = 



^2 ^^2 ( cos(/ir) + - sin(/ir)) 



(4.41) 

(4.42) 
(4.43) 



which agrees with [40H42] . Note that, as expected, the 
time and space coordinate go together forming the four- 
dimensional radius almost everywhere apart from the fi- 
dependcnt part. The corresponding zero mode densities 
also contain negative region with zeroes separated again 

by Ar = 7r/(2/i) and an increasing value at the core, 
Poo(i = 0,r = 0;//)~(l/p2 + ^2)2^ 



identifying AttTu) = v. We remind the reader that we 
started with the antiperiodic zero mode for the caloron 
and thus arrive at the antiperiodic zero mode for a time- 
dependent dyon. Therefore, its vev needs to be w; the 



transform discussed in Sec. IV is a gauge transform 
and therefore not visible in the fermionic profile, only in 
the components to be discussed now. 
The caloron's auxiliary function 



1 



r coth(47rrwT) 



COS! 



(4.45) 



is linearly divergent in d and shows the quasi-stringy be- 
haviour |35j : for large r the denominator would be zero 
on the half-line 9 = 0, if there were no exponentially small 
deviations of the coth-term from 1. The factor pv^ in 
the zero mode is 0{d) = 0{p'^). The corresponding limit 
of the Greens functions are 



lim p^f{ujT, 



T 



p—>-oo 



limpV((l 

p^OO 



lim p2/((i 

p— S-OO 



^-2TTiutT ^- lit 



X — sinh(27rwrr + ipr) 



uj)T. 



T 
'2" 



(4.46) 



(4.47) 



X — sinh(27rwrr — ipr) 



T 



1 



2 

cos(/xr) 



±i- 



(4.48) 
sin(/ir) , 



vC0sh(27rrwT) sinh(27rrwr) 



The projection matrices can be written in terms of the 
identity and the third Pauli matrix, such that the zero 
mode reads 



3. Single dyon limit p —¥ 00 

The limit of a single dyons can be obtained in the limit 
of large size p, i.e. large separation d = irp^T, of the 
dyon constituents. While r remains finite, s is large, too. 



2dr cos 9 = d — r cos 9 + 0{d ^ ) , where 9 



is the angle between r and the line connecting the dyons, 
its appearance originates from the quasi-stringy gauge of 
the caloron. 

The Greens function needed for the profile becomes 
static. 



/ 



ttT cosh(47rrTw) — cos(2^r) 
r sinh(47rrTcj) 



(4.44) 



Here, the auxiliary functions Cg and Sg have become ex- 
ponentially large in d and suppress the time-dependence 
in ip. Furthermore, Cs = Ss up to exponentially small 
corrections and these factors have canceled in all expres- 
sions. The three-dimensional Laplacian yields the profile 
from it via (4.331. It can be shown that this precisely 



equals the density Eq. (4.9 1 of the dyon zero mode upon 



T 



1 



/2nr ■\/coth(47rra;r) — cos6' 



X e 



cosh(27rra;T) 



ITS. 



sin(/ir) 
sinh(27rrwT) 

(4.49) 



Performing the time and spatial derivatives this gives 
^/TttujT 



V'o(m)^ 



2V27rr 
sm(pr) 
sinh(27rrwr) 
cos(/ir) 



Q 



i 



coth(27rrwr) 



(4.50) 

cos(/ir) 



2t:u!T cosh(27rrcjr) 
sin(/ir) 



cosh(27rrwr) 2'!TiuT sinh(27rrwT) 
T^f ■ f + tanh(27rrcjT) 



^coth(47rrtI)r) 



cost 



(4.51) 



This is the full antiperiodic zero mode of a dyon in the 
quasi-stringy gauge. The explicit appearance of and 
the angle 6 renders it different from the hedgehog gauge 
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cos(2r-/i) + cosh(27rrT)) + 



T->0 



1 / 1 I 3in^(rM)p^ ^ 



p — >■ OO 

ttT cosh(47rrTa;)— cos(2r^) 
r sinh(47rrTcij) 



- cos(27rfT) + cosh(27rrr) + sinh(27rrT) 



^ sinh(47rrra;)e' 



- cos(27rtr) + cosh(27rrT) 



1 + 



■wp-'T 



sinh(27rrT) 



- cos(27rtT)+cosh(2irrT) 



1+7^ 



COS sinh(47rrcij) + cosh(47rr(^) 

2 27rT 1 



r cos (?+coth(47rra;T) 



27rTg-Mt 



|^cos(7rtT) sinh(7rrr) cos(r/i) 



+ sin(7rtr) sin(r-^) cosh(7rr-r) 



|^cos(r/i) + |: sin(r/i) 



1 

2P 



lit ^TvitLO cos(r 
cosh(27r 



rp.) 
rTui) 



i sin(r^) 
sinh(2irrTiS) 



2p2 



cos(r^) . isin(?'/^) 
cosh{27rrTcD) ~^ siiih(2irrTcj) 



TABLE I. Table summarising some limits of functions needed to determine the zero mode density ipl{—ii)tl)o(p) = —j^dpd^f, 



and the component form given by Eq. (4.301. We labeled = f^{LjT,T/2 — ifi/2Tv 



representation used in Sec. |IV A| These only enter Q, 
and Q can be shown to be proportional to a (space- 
dependent) SU{2) matrix U 



^ ^ sinh(27rru;T) 

Q = 2 ^ ^ - U{r, 9) 

■\/sinh(47rrwT) 



(4.52) 



Performing a local (and periodic) gauge transformation 
the zero mode in hedgehog gauge finally becomes 



{e-'^'-'^^'^ {if ■ t)U'^ e"*'^^^)V'o(M) = TTwT^/^ 



cos(/ir) tanh(27rrajr) + 



27rcjr 



sin(/ir) 



(4.53) 



! ( — sin(/ir) coth(27rajr) + ^ cos(/ir) ) f ■ t 



This solution matches the one of the time-dependent 
dyon discussed at the end of Sec. |IV A| 



V. NUMERICAL RESULTS FOR THE 
STAGGERED DIRAC OPERATOR IN CALORON 
BACKGROUNDS 

In this section we will compare the caloron zero modes 
from the continuum to lattice zero modes of the staggered 
Dirac operator. The dimension of our lattices will be 
denoted by N^-Ny-Nz-Nt, the sites by an with a being the 
lattice spacing and normalized lattice vectors in the i/th 
direction by u. All observables of dimension energy will 
be measured in units of temperature T — 1/ j3 — l/{Nta), 
all length scales are given in units of /3. 



A. Discretization of the caloron incl. smearing 

The first task is to compute the lattice links Ui,{an) 
from the continuum gauge fields of calorons. To this end 



we calculate gauge transporters by discretizing the path 
ordered exponentials, 



T'exp 



A^{x)dx^ 



n 

k=l ^ 



a \ n H V 

N 



U,{an). (5.1) 



The gauge field A^{x) is singular at the caloron's center 
of mass and has large gradients on the line connecting the 
constituent dyons as one can see from Fig. 1 in |48j and 
Fig. 2 in [M] . Therefore, the number N of discretization 
points is adapted locally in steps of 40 until both the 
real and imaginary part of all matrix elements of the 
links Ui,{an) change by less than 10"^ (between N and 
3N/2). The typical value of N is 40 except at the line 
which connects the dyons. 

With this technique we obtain for a caloron with holon- 
omy parameter w = 1/4 and size parameter p = |/3 on 
a 24 • 24 • 36 • 8 lattice 184% of the continuum action 
Scant- The additional action originates from the "spa- 
tial boundary of the lattice", where onto the finite box 
we force the infinite volume calorons, whose gauge fields 
are at variance with the periodic spatial boundary condi- 
tions. To remove these artefacts we apply APE smearing 
[551 [55] with a parameter aAPE — 0.45. In Ref. [57] 
it was shown that this value provides the best match- 
ing between APE-smearing and RG-cycling. Moreover, 
in Ref. [24j APE smearing with this value has been ap- 
plied to SU(2) Monte-Carlo gauge configurations to re- 
veal their dyon content. 

After 275 APE smearing steps the boundary arte- 
facts are sufficiently smoothed out reaching now 1.07% of 
Scant- In order to have comparable results all gauge links 
are smeared until the lattice action reaches 1.07 Scant- 
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B. Staggered Dirac operator at finite cliemical 
potential and its symmetries 

We use the standarcQstaggered operator, which in our 
notation has the form 

Difi; n\m) = ^J2 vAn){Sn+au,mUAan)e''^^-'^^ (5.2) 

U — l 

-5n-a^,™C/^(am)e-'^^^-^}. 

where rii,{n) = (— i)"i+ - "'^-i is the staggered sign. 
At vanishing chemical potential D is anti-hermitian 
with purely imaginary eigenvalues as in the continuum, 
whereas at nonzero chemical potential the hermiticity re- 
lation Eq. (2.4) to opposite ^ holds. We have imple- 



mented the chemical potential by exponential factors on 
all temporal links. 

This operator obeys a remnant of chiral symmetry 



{D{fi),7j5}^0 ,75(«)-(-l) 



(5.3) 



yielding eigenvalues ±A as in the continuum. The anti- 
unitary symmetry, however, differs from the one of the 
continuum Dirac operatoij^ It reads 



D*{^i) = TlDi^i)T2 



(5.4) 



with T2T2 = —1. The staggered Dirac operator is thus 
in the universality class of the Gaussian symplectic en- 
semble (the expected transition to the continuum sym- 
metry has been investigated in [SSI [SS]). As a conse- 
quence, T2'il^* is another eigenmode wit h eig envalue A*. 
Together with the chiral symmetry, Eq. (5.3), this yields 
(Kramer's) degenerate eigenvalues for vanishing chemical 
potential. At nonzero chemical potential the eigenval- 
ues come in quartets {A, —A*, —A, A*} with eigenmodes 
{V", %T-2 ■'/'*, %V': T21P*}, respectively. 

On our smooth backgrounds an approximate four-fold 
degeneracy is expected from the four tastes of the stag- 
gered operator. Thus, instead of each exact continuum 
zero mode a low-lying quasi-zero quartet shall occur. At 
^ = topological quasi-zero modes of the staggered op- 
erator are known to exist. The chirality is represented 
by the four- link operator Esq, with which these modes 
have expectation values close to ±1, in contrast to higher 
modes |60H63) . Both, nonzero real parts and nonzero 
imaginary parts of the eigenvalues could be induced by 
discretization and finite volume effects. They may also 
differ for quasi-zero eigenvalues and singular values. 

On the smeared configurations we find modes separated 
by orders of magnitude from the rest of the spectrum (see 
Eig. |6|. We investigate the profiles of these quasi-zero 



Improvements are not necessary on our smooth backgrounds. 
* because no charge conjugation matrix is needed for the staggered 
operator being a scalar in spin spcae 



modes without further analysis of their continuum and 
infinite volume limit. As a side remark we state that 
on a 48 • 48 • 72 • 16 lattice the smallest eigenvalue of 
a caloron with maximal nontrivial holonomy and zero 
chemical potential reach zero within double precision. 

As the staggered operator and its eigenmodes cannot 
be made real as in the continuum, the argument of the 
reality of the zero mode profile needs to be revisited. 
Erom the computations it will turn out that the stag- 
gered quasi-zero modes have real parts vanishing to ma- 
chine precision, such that the quartets practically con- 
sist of two degenerate pairs, (A,— A*) and (—A, A*), on 
the imaginary axis. If an eigenvalue A is degenerate with 
—A* the corresponding subspace is spanned by 

(4>{x;fi),r]r,T2ip*{x;^)j = (^^i(a;; ^), ■02(a;; m)) (5.5) 

In Appendix [B] we show that a particular linear combi- 
nation of these modes obeys the bi-orthonormalization 



d''x^lj{x;-fiyjN{x;fi) = Smn, M,Ne {1,2} (5.6) 



and that the profile averaged over these two modes 
I 2 

Poo{x;n) = 2 ^ ^M(^;-M)V'M(a;;Ai) (5-7) 

Af=l 

is real. Consequently, we will compare the real profiles 
(5.7) to those in the continuum. 



C. Results for eigenvalue spectra 



On the configurations described in Sec. V A we mea- 
sured 256 eigenmodes with smallest magnitude by virtue 
of ARPACK t64j, Armadillo [65. and the G++ Boost Li- 
braries. As expected the eigenvalues of the staggered op- 
erator become complex when switching on the chemical 
potential, shown in Eig. [6j The quasi-zero eigenvalues, 
however, remain basically unchanged at finite chemical 
potential. Erom that figure one further recognizes, that 
the other low-lying eigenvalues of the caloron are similar 
to generalized Matsubara frequencies. By that we mean 
the Dirac spectrum in free backgrounds with constant 
Polyakov loops of the same holonomy as the caloron 



Ut = exp (27ri diag(w, —uj)/Nt) 



(5.8) 



and trivial space-like links, = Uy = = 1, such that 
the action vanishes. This background gives rise to the 
following staggered eigenvalues 



(5.9) 



S {w. 

i—x,y,z 



2tt 



{kt±uj 



if — ifJ^/T 
2^ 
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FIG. 6. Spectrum of the staggered Dirac operator in a caloron 
background of maximally nontrivial holonomy (oj = 0.25, p = 
0.75/3) and different chemical potentials on a 24 ■ 24 ■ 36 • 
8 lattice. Note that in order to be able to compute more 
eigenvalues, this lattice is smaller than the one used for the 
profiles later in Fig. |8] For comparison the triangles show the 
eigenvalues of a free Polyakov loop background with maximal 
nontrivial holonomy (see text). These are highly degenerate. 



with fc^ e (— . . . , the wave numbers of the 



plane wave eigenmod.es, for vanishing fi cmp. 



For 



low-lying modes, ^ N^, the sines disappear and the 
eigenvalues follow the continuum dispersion relation. 

The similarity of the caloron spectra to these general- 
ized Matsubara frequencie^ holds inch nonzero chemi- 
cal potential, cf. Fig. |6] The caloron's quasi-zero modes, 
which are of topological origin, are of course not reflected 
in the Matsubara frequencies. 

An interesting effect occurs in the low-lying spectrum 
when the boundary condition phase ip equals the holon- 
omy w, such that in the free case they compensate in 
one color sector (this is the case where the caloron's 
zero mode cannot be analytically continued since at this 
boundary condition it hops between the constituents) . At 
vanishing chemical potential the free spectra then possess 
a zero mode for k = = kt- With nonzero /i, the eigen- 
value with these momenta becomes 

T 



± Nt sinh 



Nt J 



± 



(5.10) 



and hence is located on the real axis. With the lowest 
frequency kt = 0, and nonzero fc, there is a tower of 
states, which at /i = form the low-lying modes on the 
imaginary axis, while at nonzero /i they become 



A(/i) 
T 



= ±Nt, sinh-" 



Nt 



sm 



2n 

N 



h] (5.11) 



^ A similar similarity holds for the eigenvalues of the gauge covari- 
ant Laplacian |66| . 



FIG. 7. Spectrum of the staggered Dirac operator in a caloron 
background of maximally nontrivial holonomy (same as in 
Fig. [6| and intermediate boundary condition (periodicity up 
to i) plus the corresponding free spectrum. The lowest eigen- 
values are either real or purely imaginary, see text. 



In other words, these modes stay on the real axis with 
magnitudes smaller than the one above and then become 
purely imaginary with growing magnitudes. How many 
of these either real or purely imaginary eigenvalues oc- 
cur on the lattice depends on the aspect ratios Ni/Nt- 
Only when moving to the next frequency \kt\ — 1 or to 
the color sector, where uj does not compensate but adds 
to the boundary condition, do these eigenvalues become 
fully complex. 

The lowest staggered eigenvalues of the caloron in this 
case are also either real or purely imaginary, see Fig. [7] 
for maximally nontrivial holonomy and boundary phase 
if — IT /2 (i.e. periodicity up to a factor i). One of the 
caloron's real eigenvalues is again close to its free coun- 
terpart (5.10). 



D. 



Comparison with the exact zero mode for the 
caloron 



Having found quasi-zero modes in the staggered spec- 
trum, we expect the corresponding eigenmodes to asymp- 
tote to the continuum zero modes of Sec. |IV C[ To get 
the eigenmodes to be bi-orthonormal, we have computed 
eigenmodes at /i and —fj, and used the linear combina- 
tions defined in Eq. (5.7). The comparison is shown in 
Fig. [8] for fixed chemical potential, /i = 2AT, while the 
holonomy increases from trivial to maximally nontrivial. 
To compare the numerical solution to the continuum one, 
we fix the time t and the coordinates z (along the caloron 
axis) and y (perpendicular to the caloron axis) as close 
as possible to the twisted dyon's core, on which the (an- 
tiperiodic) zero mode is localized. Our lattices are shifted 
such, that the caloron axis runs through the middle of a 
plaquette. Therefore, we choose y = t = a/2 = /3/{2Nt). 
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In order to show two-dimensional plots, we finally fix the 
z-coordinate. We choose that value of z, for which the 
numerical density poq{x; ^jl) (cf. Eq. (5.7)) as a function 
of the remaining coordinate x is maximal. The analytic 



density (4.33 4.34[ )) is plotted along the corresponding 
continuum line. For holonomics cj = 0.0, 0.125, 0.25 this 
amounts to shifts in z of 0.5a, 0.8a and 1.1a from the 
formal dyon core r — Q towards the other dyon. 

We find that the numerical densities coincide well with 
the continuum zero mode densities. In particular, the 
negative regions are confirmed by the numerical solutions 
(cj = 0.125, 0.25 in Fig. [s]). That for increasing holon- 
omy the mismatch becomes stronger can be explained by 
considering the action densities of those calorons. For 
Lj — Q there is just one dyon which is far away from 
the boundaries of the lattice, while in the case of maxi- 
mal nontrivial holonomy two separated dyons are in the 
same box. This means that the finite volume effects are 
stronger in the case of maximally nontrivial holonomy. 



VI. OVERLAP MATRIX ELEMENTS 



A very important quantity for building a dyon liquid 
model in analogy to the instanton liquid model is the 
overlap matrix element [HI |42] . 



I. J 



d'^x tp\{x; + ifi'-fo)U'il;j(x; fi) , (6.1) 



where ipij^j are zero modes of dyons labeled by the 
index / and of antidyons labeled by the index J, re- 
spectively. This form reflects the residual U{1) gauge 
freedom, as the two dyons can always be superposed in 
gauges differing by a quasi-abelian gauge transformation 
U = expiiatb ■ f) which preserves the holonomy. This 
replaces the general SU{2) relative gauge transformation 
in the case of the instanton anti-instanton overlap matrix 
element (see ^^)- 



Our formulas Eq. (4.301 apply only for the case that 
holonomy is chosen along the line that connects two con- 
stituent dyons, and that this line is chosen along the z- 
direction. This also applies to the single dyon limit of 
the previous subsection. It is easy to generalize this ex- 
pression to the arbitrary direction of the Dirac string 



FIG. 8. Comparison of the analytical and numerical density 
of the caloron zero mode at ^ = 2.4 T, on a line perpen- 
dicular to the caloron axis and intersecting with it near the 
twisted dyon (for more details on the chosen coordinates see 
text). All panels correspond to caloron with size parame- 
ter p — 0.75/3. The holonomy increases from top to bottom: 
uj — 0.0, 0.125, 0.25. The black crosses are from numerical 
solutions on a 36 ■ 36 • 54 • 12 lattice, while the red solid line 



correspond to the analytical solution in Eq. (4.331 



■0/(t,r;/i) 



T 



2V27rr ^coth(47rrwT) - cos 6*1 
cos(r/i) 



ISl ■ T 



cosh(27rra;r) 
sin(r/i) 



sinh(27rrti;r) 



^-^^^a,e} , (6.2) 
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FIG. 9. Plots of numerical integration of functions T^j , Eq. (6.71, as a function of distance d in units of v. The plots show 
results for jj, — (solid), fi = v (dashed) and n = 2v (dot-dashed) . Note that the behavior is similar to that of the behavior of 
two spatially separated instantons (see Fig. 1 in |14| ) 



2\/27rr ^coth(47r swT) - cos 02 



+ tS2 ■ T — 



vCOsh(27rsa;r) 
siii(s/j,) 



sinh(27rswr) 



where now 9i^2 are the angles between r 
^ Si_2 respectively (see Fig 



string^l 

superpose in the algebraic gauge in which — >■ away 



a,e} . (6.3) 



s and the Dirac 
lOl). Note that we 



from the dyons, cf. Sec. IV B 1 The dyons in question 



are those that carry the anti-periodic zero modes (the 
"heavy" ones), although the same procedure can easily 
be adapted to the case of dyons carrying the periodic zero 
modes. 

However, the above solutions assume holonomies in the 
directions si and —S2 for the dyon and antidyon. To 



S2 



FIG. 10. Configuration of dyon (shaded) and antidyon 
(white) and their Dirac strings (dashed lines) discussed in 
the text. 




Note that because for a selfdual field A^(t,x), we have that 
Af^{t,x) = (Ao(t, — x), — A(t, — a;)) is anti-selfdual. That means 
that if Ip{A)ip = then ]p(A)'yoip{t, —x) = 0, so one can con- 
struct an antidyon (anti-caloron) zero mode solution from the 
dyon (caloron) zero mode solution by flipping x — > —x and mul- 
tiplying by 70. Note that when we construct the antidyon zero 
mode in this way we also have to flip the direction of the Dirac 
string. 



compensate that we define rotation matrices 
UiSi ■ tU\ ~ lj ■ t , 



U2S2 ■ tU, 



(6.4) 
(6.5) 



where uj is an arbitrary unit vector parametrizing the fi- 
nal direction of the Polyakov loop in color space (which 
without loss of generality can be chosen in the third di- 
rection). So what we need to compute for the overlap 
matrix element is 

Tjj— / d^xipl{x;—fi)l{lui{i{~i0 + ifijo)ipj{x;fi) 



(6.6) 

where U = exp^iaio-f) is the relative quasi-abelian gauge 
transformation between the dyons mentioned already at 
the beginning of this section. Note that for instanton a 
relative SU{2) gauge transformation was relevant. Here 
the relative quasi-abelian gauge transformation is rel- 
evant, while the rest of the subgroup SU{2)/U{1) re- 
duces to the orientation of the strings (i.e. angles between 
strings and holonomy). 

We leave the problem of the hopping element of a 
caloron in full generality for future work, but here to illus- 
trate the qualitative behaviour let us take a simple con- 
figuration of the dyon-antidyon system with their strings 
lying on the z-axis, and pointing away from each other. 

Writing the overlap matrix element 

Tjj^i {T]j cos a -f TIj sin a) (6.7) 

it is not very difficult to see that T^j = 0, at /i = 0. 
This happens simply because of the index structure of 
the above expressions, resulting, upon index contraction, 
in vanishing contribution. For /i 7^ this is no longer the 
case, as there is an additional term propo rtiona l to T3 in 
the expression for the zero mode (see Eq. (4.49)). In Fig. 
|9] and [Tl] we show the results of numerical integration 
of the overlap matrix element as a function of distance 
between the dyon and anti-dyon. 
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FIG. 11. Logarithmic plots of T^'? as a function of distance d in units of v, for = (solid) 
(dot-dashed). The dashed, straight lines are linear functions with slope —1/2 to guide the eye. 
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V (dashed) and fj, = 2v 



Because of exponential localizations, the most domi- 
nant contribution to the overlap matrix at large distance 
would naively come from the two regions where the zero 
modes are localised. However since they decay exponen- 
tially, along the line connecting the dyon and antidyon, 
there is a significant contribution in this region, turning 
the asymptotic behavior into exponential ^ exp(— u(i/2) 
(with no algebraic factors). 



VII. SUMMARY 

We have analyzed zero modes at nonzero chemical po- 
tential in topological backgrounds. Starting with non- 
trivial holonomy calorons, for which we have analytically 
continued known zero modes from imaginary chemical 
potential (i.e. phase boundary conditions in the tempo- 
ral direction), we have derived the zero modes for trivial 
holonomy calorons and instantons known in the litera- 
ture as well as for dyons, for which a spherical ansatz 
was explored, too. All of these procedures have a direct 
generalization from SU(2) considered here to SU(N). Zero 
modes present at /i = due to index theorems remain at 
nonzero /i. 

The most prominent features of these zero modes are 
higher values at their centers and negative regions in their 
(pseudo-) density. In dyon backgrounds the former in- 
creases asymptotically like /Lt^, while the distance between 
the zeros between regions of different sign is proportional 
to l//i. The envelop of the decay is to leading order in- 
dependent of /i. 

On lattice-discretized calorons we have found com- 
pletely analogous quasi-zero modes for the staggered 
Dirac operator. For the latter this means that although 
exact chiral symmetry is lacking, the topological part of 
the continuous spectrum on such smooth backgrounds 
is well-reflected (in eigenmode quartets) also at nonzero 
chemical potentials. Therefore, this numerically cheap 
operator could be sufficient for the analysis of ensembles 



of topological ensembles (or realistic lattice Monte Carlo 
configurations) at nonzero density. 

The overlap matrix elements of such zero modes are a 
key ingredient to the interaction of such topological ob- 
jects. As we have demonstrated by virtue of numerical 
integration, the exponential decay of the overlap matrix 
elements with the distance are the same as for zero /i, 
but their alternating behavior will cause an additional 
wash-out effect similar to that of instantons [T3]. We 
also showed that a term in the overlap matrix element 
appears, that is related to the relative abelian orienta- 
tion of the dyons and that vanished at /i = 0. More 
consequences of the overlap matrix elements at nonzero 
chemical potential are left for future studies. 
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Appendix A: Reality of profiles in the continuum 

If 0(/i) vanishes on the zero mode -00 (/i), then from 
the co njuga te equation and the anti-unitary SU(2) sym- 
metry (2.5) it follows that 0(/i) vanishes on WipQ{fj,)* as 
well. For a nondegenerate zero mode the latter must be 
proportional to the zero mode up to a factor 



(Al) 



This factor a is indeed a phase, which can be seen by 
acting with W on the complex conjugate equation using 
WW* = 1 and the previous equation again: 



ww*M^^) 

■00 (m) 



Wai^lrM^^r 

a(A*)*a(M)V'o(M) 



(A2) 
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Next we prove that a{—fi) = a(fi). For that we use the 
bi-orthonormahzation at ±fi, — W and (Al) at — /i 
to write 

a{fj.) = / d'^xil)l{^i)a{-^i)ipQ{-^j) = a{-^j) (A3) 



Finally, from (Al), its hermitian conjugate and the uni- 
tarity of W: 



'05(-m)V'o(m) 



(A4) 



which just means that the density pooifi) is real for all /j, 
and X. 



Appendix B: Bi-orthonormalization and reality for 
degenerate staggered eigenmodes 



With the definition (5.5) the densities in the bi- 
orthonormalization 



V'L(a;;-Ai)'0w(a;;Ai) = Pmn{x]^i) (B1) 

(B2) 



obey 

and have the form of a quaternion 
/ 



p(x;ij) 



jd^^pi^-.p) 



g(x]p) h{x;p) 
\ -h*{x;p) g*{x;p) 

G{p) Hip) 
H*{p) G*{p) 



(B3) 



(B4) 



i.e. the modes are not necessarily bi-orthogonal (unless 
H = 0). 

Choosing two different modes in the subspace amounts 
to a right multiplication 

i^M{x\p)^i^N{x-p)CNM{p) C{p)eGL{2,C) (B5) 

where the coefhcient matrix C{p) may depend on p, but 
must be a constant wrt. x. The density transforms as 



p{x;p) C'^{~p)p{x]p)C{p) 



(B6) 



preserving the property (B2). The integral transforms 
accordingly, 

fd'^xp{x]p)~^C\-p) jd^xp{x]p)C{p) (B7) 
With the choice 

C{p) = \ fd'^xp{x;p)] ^'"^ = G\-p) (B8) 



(B9) 



one achieves the desired bi-orthonormalization 



d^x p{x] p) — > Ij 



and for the sum of the densities 



^i'\i{x]-p)'4)M{x;p) -> irC'^{-p)p{x]p)C{p) 



M=l 



= tr 



a H 



-H* G 
Gg* +Hh 




|G|2 + |i?|2 



-I- c.c. 



(BIO) 



which is real. That this expression is the sum of two 
densities is accounted for by a factor 1/2 in Eq. (5.7). 
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